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Abstract. We study the critical behavior of the Anderson localization-delocalization transition in corner¬ 
sharing tetrahedral lattices. We compare our results obtained by three different numerical methods namely 
the multifractal analysis, the Green resolvent method, and the energy-level statistics which yield the sin¬ 
gularity strength, the decay length of the wave functions, and the (integrated) energy-level distribution, 
respectively. From these measures a finite-size scaling approach allows us to determine the critical param¬ 
eters simultaneously. With particular emphasis we calculate the propagation of the statistical errors by a 
Monte-Garlo method. We find a high agreement between the results of all methods and we can estimate 
the highest critical disorder Wc = 14.474(8) at energy Ec = —4.0 and the critical exponent u = 1.565(11). 

Our results agree with a previous study by Fazileh et al but improve accuracy significantly. 

PACS. 71.30.+h Metal-insulator transitions and other electronic transitions - 72.15.Rn Localization effects 
(Anderson or weak localization) - 64.60.F- Equilibrium properties near critical points, critical exponents 


1 Introduction 

On a very fundamental level the Anderson model of local¬ 
ization [1] demonstrates the influence of quantum interfer¬ 
ence on the diffusion of an electron within a lattice. The 
variation of disorder, e.g. by a specific choice of random 
on-site potential energies, leads to a phase transition of 
the electronic wave functions between localization and de¬ 
localization (LD) corresponding to insulating and metal¬ 
lic behavior, respectively. The critical behavior at the LD 
transition is described quantitatively by the scaling theory 
of localization US]. In particular, the localization length 
and the conductance behave according to a power law in 
the vicinity of the transition. The corresponding critical 
exponents are expected to adopt universal values as pre¬ 
dicted by random matrix theory m meaning their values 
depend only on the classification according to the symme¬ 
try of the underlying Hamiltonian. Universal results ob¬ 
tained for the rather abstract Anderson model of localiza¬ 
tion are therefore applicable also to complex experimental 
measurements [6]. 

The Anderson LD transition was studied extensively 
during the last decades (for reviews see [7l[8l[9]). A major 
subject of theoretical research covered the numerical ver¬ 
ification of the concept of universality. Here considerable 
progress was achieved during the last decade due to high 
accuracy calculations combined with extensive data and 
error analysis maim- 

In this work we study the Anderson LD transition in 
corner-sharing tetrahedral (CST) lattices. The CST lat¬ 
tice appears in normal spinel structures as a non-bipartite 


sublattice, e.g. the Ti atoms in a LiTi 204 crystal. Gradual 
replacement of Ti atoms by Al atoms leads to a transi¬ 
tion from metallic to insulating behavior. Recent numer¬ 
ical analyses mm of this LD transition focused on an 
estimate for the critical concentration of Al. The results 
were inconsistent with experiments which was attributed 
to the use of one-electron models missing electron-electron 
correlations mm- Nevertheless the LD transition in a 
CST lattice should obey the predictions of random matrix 
theory. According to the symmetry of the corresponding 
Anderson Hamiltonian the model belongs to the Gaussian 
orthogonal ensemble (GOE). For GOE in three dimen¬ 
sions the most accurate estimate of the critical exponent 
u of the localization length m is obtained for a simple 
cubic (sc) lattice by the transfer-matrix method [15] . The 
result of u = 1.571 ± 0.008 agrees well with other recent 
results from multifractal analyses of the wave functions 
[IOllIl|T7|. Since the symmetry of the Hamiltonian is in¬ 
dependent of the crystal lattice itself, the same universal 
exponent u is expected also for other lattice geometries, 
which was demonstrated for the face-centered-cubic (fee) 
and the body-centered-cubic (bee) lattice [18] . In this pa¬ 
per we investigate the critical behavior, i.e. the exponent 
u of the localization length, for the LD transition in a 
CST lattice. We compare results of three different meth¬ 
ods. Primarily we use the multifractal analysis (MFA) to 
obtain the phase diagram. This method utilizes multi¬ 
fractal fluctuations of the wave functions at the critical 
point. Secondarily we apply the Green-resolvent method 
(GRM) to calculate the exponential decay of a wave func¬ 
tion along a bar-shape geometry. The commonly preferred 
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transfer-matrix method (TMM) is not applicable in our 
case because of the connect ivit}Q of the CST lattice. As a 
third independent approach we study the distribution of 
energy-level spacings (ELS). Because of the finite system 
size of the calculations all three methods require a finite- 
size-scaling (ESS) procedure in order to obtain the critical 
parameters. 

This paper is organized as follows. In Sec.[^ we shortly 
introduce the Anderson model of localization and describe 
the Hamiltonian for the CST lattice. In the following Sec¬ 
tion we explain the MEA, GRM, ELS, and the ESS ap¬ 
proach. In Sec. we present our numerical results and 
compare the methods. Einally, in Sec. we conclude with 
a short summary. 


2 Anderson model of localization and the 
corner-sharing tetrahedral lattices 

We consider the Anderson Hamiltonian for non-interacting 
electrons in site representation 

H = Y^Vi\i) {i\-^\i) {j\ . (1) 

Vi are randomly chosen on-site potentials distributed uni¬ 
formly in the interval [—IT/2, IT/2]. The interval width 
IT parametrizes the disorder strength. Every site is con¬ 
nected to its six nearest neighbor sites according to the 
structure of the CST lattice as shown in Eig. 

The CST lattice is built from regular tetrahedrons. 
The shared corners represent the lattice sites and the edges 
the nearest neighbor connections. The different colors dis¬ 
tinguish tetrahedrons that can be transformed into each 
other by simple translations. The CST lattice is a mem¬ 
ber of the fee Bravais-lattice class. We build up the lattice 
using primitive unit cells (UCs) as the smallest possible 
building block as illustrated in Eig.[^ It contains four lat¬ 
tice sites. The six edges of the blue tetrahedron represent 
the intra UC connections. Taking into account the trans¬ 
lational symmetry, there are also six different inter UC 
connections visualized by the edges of the green tetrahe¬ 
drons. The three Bravais-lattice vectors span a rhombohe- 
dron. Because the CST lattice is a sublattice only, half of 
the rhombohedron is empty. Note that alternatively the 
CST lattice can be constructed also in a non-primitive 
way from a fee lattice. Such a UC contains 16 sites. In the 
following we will express the system size always in units 
of primitive UCs stacked along the three Bravais-lattice 
vectors. In general we use periodic boundary conditions. 


^ For all crystal orientations the construction leads to singu¬ 
lar link matrices between successive layers. Also, a subdivision 
of a layer into several sublayers as it was done in triangular 
lattices m would not help. 



Fig. 1. A part of the CST lattice. The spheres embody the 
lattice sites and the edges of the tetrahedrons epitomize the 
nearest neighbor connections. The coloring of the tetrahedrons 
indicates the two translation invariant groups of tetrahedrons. 



Fig. 2. Primitive UC of the CST lattice. The shaded rhombo¬ 
hedron represents the Bravais lattice, the red spheres the basis 
sites. The gray spheres are sites of neighboring UCs. On each 
corner of the rhombohedron is one green tetrahedron (not all 
shown). Only the visualized green tetrahedrons have a connec¬ 
tion with the blue tetrahedron in this UC. 

3 Methods 

3.1 Finite-size-scaling approach 

The LD transition is a second-order phase transition char¬ 
acterized by the algebraic divergence of the correlation 
length oc \uj\~^ as the transition is approached {uj 
0). The dimensionless scaling variable uj = defines 

the distance from the critical point pc? in onr case in 
terms of either energy (p = E) or disorder strength (p = 
IT). This behavior is observed only for an infinite system 
while our calculations are limited to finite system sizes. 




Martin Puschmann et al: LD transition in CST lattices 


3 


We therefore employ a FSS ansatz in order to extract the 
critical parameters of the transition. The implementation 
is based on an expansion function for the measured quan¬ 
tity close to the critical point and considers both relevant 
and irrelevant influences [20]. Irrelevant contributions are 
quantified by an irrelevant exponent y describing finite- 
size corrections that eventually disappear with increasing 
system size. Let us consider e.g. the localization length 
A{uj^L) for a system with characteristic length L. Nonlin¬ 
earities are taken into account by expanding the uj depen¬ 
dence in terms of Taylor polynomials and the dependence 
on L in terms of Chebyshev polynomials of the first kind 
For the expansion orders mi, and m[ the 
scaling function is given by 


A = F 

with 


f2,Li,f2iL-y = Y2 - Pr,. 


n=0 


n\ 




( 2 ) 


QrL-^ 


^ ^ Q'n,feF'/e (i2r-L ^ ) , 

k=0 

( 3 ) 

CJ + y , and 

. n! 

n=2 

mi 

( 4 ) 

^ n\ 

n=l 

( 5 ) 


is given by the sum over the squares of the corresponding 
wave amplitudes of the eigenstate </> 


iGbox h 


The mass exponent for moment q is defined by 


rg{(p,l,L) = |im 


lnPg{cP,l,L) 
In A 


(6) 


( 7 ) 


with the generalized mass 


b 

and relative box size A = l/L. According to Eq. 0 the 
value of Tq can be extracted as slope from the linear re¬ 
gression of InP versus In A with emphasis on small box 
sizes. In this work we compare several ranges between 
0.1 < A < 0.5. 

The singularity spectrum f{o:q) [23| contains further 
information about the wave function and can be calculated 
by a parametric representation of the singularity strength 


dq 


'Fl = lim 
dq A^o In A 


and the fractal dimension 


( 9 ) 


This expansion contains = (ni -h l)(nr +1) Tmi + 
1 + 0(ni) free parameters to be determined by a nonlinear 
fitting procedure. 

In order to estimate the errors associated with the non¬ 
linear fit reliably, we use a Monte-Carlo method [10]. We 
generate at least 10^ synthetic datasets from our original 
data by adding Gaussian noise to each data point. The 
individual error of each original data point defines the 
standard deviation of the noise. Combining the fit results 
for all synthetic datasets we obtain distributions for the 
critical parameters. In general these distributions have a 
Gaussian shape which allows us to identify standard de¬ 
viations easily. In contrast to this stochastic inaccuracy it 
is more complicated to identify systematic deviations. We 
use several expansion orders to reduce the influences of a 
specific fitting function. In some cases we observe devia¬ 
tions from Gaussian distributions, which we attribute to 
local minima of the sum in the nonlinear fit. Their oc¬ 
currence can be reduced by a reasonable choice of smaller 
expansion orders. 


3.2 Multifractal analysis 

In the vicinity of the LD transition multifractal fluctua¬ 
tions of the wave functions are observed [HEa.The MFA 
quantifies these fluctuations by a box-counting algorithm 
using different moments of the probability distribution. 

A system of UCs is subdivided into boxes of size . 
The probability to find an electron inside a specific box 


with 


fq — Q.(^q 


lim 

A^O 


In A 


-q ^ l^b 


( 10 ) 


( 11 ) 


Numerically fq and aq are determined in the same way as 

In the standard box-counting algorithm the system is 
divided into boxes without overlap. Thus only box sizes 
with an integer ratio for L/l are allowed. Depending on 
L the number of possible A values can be very restricted. 
This limitation can be eliminated by using periodic bound¬ 
ary conditions together with a folding-back of protruding 
sites into the original system mEM- 

A finite single wave function is insufficient for a sta¬ 
tistical statement. Therefore we use a set of wave func¬ 
tions and averaged quantities. Two different averages are 
commonly used, the ensemble (ens) and typical (typ) av¬ 
erage [Ol fTbllTT] . Exemplified for the mass exponent one 
obtains 


and 


ens 

‘q 


lim 

A—^0 


In A 


T^P = lim 
® A-!-0 


In A 


( 12 ) 


(13) 
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with arithmetic average (•)0 over all considered wave func¬ 
tions. The averaged quantities from Eqs. § and ( [Tq| ) can 
be derived analogously. 

For application of the MFA to the CST lattice, our sys¬ 
tem is stacked with primitive UCs and thus has a rhom- 
bohedral form. Therefore it is convenient to use rhombo- 
hedrons as boxes. The shape has no further effect on the 
algorithm. The primitive UC is used as the smallest possi¬ 
ble box. A further subdivision seemed to be inappropriate, 
because empty boxes would follow. 


3.3 Green resolvent method 


The GRM is based on a iterative formulation of the time- 
independent Green function G = (Z — for a system 

described by a Hamiltonian H [25ll26[l27| . Z = (F + ±r])l 
is a diagonal energy matrix with the complex quantity 77 . 
The system has a bar-shaped geometry with cross section 
L X L and length N L and grows layer by layer at one 
end of the bar. This process is considered as a perturbation 
A of the previous system H and thus expressed with the 
Dyson equation 

G = G + GAG . (14) 


G is the corresponding Green function of H. The further 
description uses a notation where G • ^ ^ is the submatrix 
of a system with length N which couples the ith and jth 
layers. In the same manner the Hamiltonian can 

be represented as tridiagonal block matrix with the layer 
Hamiltonians as diagonal entries and the inter¬ 
layer connections or as nondiagonal entries. 

The matrix without the last row and column is the Hamil¬ 
tonian of the previous step. The last nondiagonal 

block matrices V n,n-\-i and V n-\-i,n are used for the per¬ 
turbation A of the matrix is the com¬ 
bination of the matrix and the last diagonal entry 

^TV+i^AT+i* ^ow Eq. Hdk can be simplified to the four 


cases 



'^Ar+i,Ar+i 

— '-^Ar+i,Ar+i y n+i,n^n,n+iJ 

(15) 

^N^l,k<N 

— ^N+1,N+1^ N-\-l,N^N,k ’ 

(16) 

^j<N,N+l 

= G^/^Vn,n+iG^n+i,n+i ’ and 

(17) 

j<N,k<N 

- '^j,k N,N+l'^M+l,k ■ 

(18) 


Combining Eqs. ([l^ and (17) then yields 


G^n+i]n+i = (Z - Hjv+i,iv+i - Siv)“' (19) 


with Tin =Y iv+i,ArG^]^VAr,Ar+i- For the analysis of the 
critical behavior we employ the localization length A of 
the wave function for a long bar-shaped geometry 


1 

A 


lim lim Intr 

2N rj^O N^oo 




— — lim lim In 

N 77^0 N^oo 




( 20 ) 


with Frobenius norm ll'Hp. The submatrix 


gS = Gd li ^M+iG^yd+i (21) 

i=l 

can be calculated iteratively by using Eq. 0 - In order to 
stabilize this calculation numerically MacKinnon et al. m 
used a transformation on simple cubic lattices, which can 
be generalized by the following steps where the size k re¬ 
flects the number of layers comprised in one step: 
i) The calculation is initialized by 




r-(i) 

^ 1,1 


k-l 


n 






i=i 


( 22 ) 


with bi = ||Bi||p. Bi = ^ is the normalized matrix of 
Bi. ii) At each further step we build 

k 

B,+i = B, n l^o+iGkti+i (23) 

i=i 

and Bi+i = with 6^+1 = ||Bi+i||p. The norms bi 

are related to the corresponding Green functions for an 
integer number of steps N/k by 

(24) 


lim lim (In b) 

? 7^0 N^oo 

(25) 

Here (•) describes the average over N/k steps. Correspond¬ 
ingly, the statistical error of A follows from 

A^ 

var A = var ((In 6 )) . (26) 

rv 

The step size k is limited by the two following facts. On 
the one hand it should be small {k N) with regard to 
the numerical stability. On the other hand it should be 
large (/c ^ 1) to get an accurate error estimation. This is 
due to the fact that successive bi are correlated to each 
other m- That correlation will be reduced by increasing 
k. We use k = 1000. 




N/k 

t[bi 

i=l 


Eq. (20) then transforms to 

11 N/k ^ 

— = — — lim lim In 6 ^ = — — 
A Nv^ON^oo^^ k 


-1 


3.4 Energy-level spacing 

In contrast to the previous methods the ELS approach 
relies on the eigenenergies instead of the eigenstates of 
the system. We consider a subspectrum, which consists 
of the + 1 adjacent eigenvalues Ck closest to energy E. 
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Fig. 3. Analytical level-spacing distributions for the extremely 
localized (Poisson distribution) and the extremely delocalized 
case (Wigner surmise) in comparison to numerical results for 
three disorder strengths near the transition with energy Ec = 
—4.0. The system size is L = 50. 

The subspectrum is then unfolded by the normalized level 
spacings to obtain a dimensionless quantity 

For the ELS distribution p{s) = {6{s — s'))^, it follows 
that (1)^ = 1 and {s)^ = 1. The expected distributions are 
shown in Fig.|^ For the extremely localized case the eigen¬ 
values are uncorrelated and the ELS distributions follow 
the Poisson law = exp(—s). In the delocalized case 

levels repel each other and the spacings are distributed ac¬ 
cording to the Wigner surmise p{s)^^^ = exp( —Js^). 
All ELS distributions cross in a common point near sq = 
2.0 which is utilized in the definition of a new measure 

m 

So 

■j = J dsp{s) = {0 {sq - s'))^, . ( 28 ) 

0 

The average value of 7 over various realizations of the sys¬ 
tem for a given disorder strength W allows us to increase 
the accuracy of our result and to estimate its statistical 
error by the standard deviation. 

4 Results 

First we discuss the density of states (DOS). The DOS 
is calculated by diagonalization of the Anderson Hamilto¬ 
nian 0 . Its dependence on W is shown in Fig.[^ Although 
the number of 6 nearest neighbors is the same as for the 
sc lattice, the DOS is very different. For systems without 
disorder (W = 0 ) we obtain four separate energy bands 
due to the four basis sites. Two flat bands are degener¬ 
ate and engender the delta peak at energy E = 2.0, i.e. 
at the high energy band edge. The other two bands are 
distributed symmetrically with respect to E = —2.0. The 
total DOS spans the interval E G [—6, 2] and is asym¬ 
metric in shape, which follows from the non-bipartiteness 



Fig. 4. DOS for a system of 19 x 18 x 17 UCs for different disor¬ 
der strengths W. The results are averaged over 200 realizations 
of disorder. The bin width is 0.1. 


of the CST lattice. The energy range ensues from the fee 
Bravais-lattice class and the number of nearest neighbors. 
With increasing disorder strength W the bands broaden 
as can be seen in Fig. For W = 2 the two dispersive 
lower bands are not very much affected, their DOS struc¬ 
tures are only slightly washed out, while the two flat bands 
are strongly influenced and the delta peak is significantly 
widened to the interval E G [1,3] corresponding to the 
width of IF = 2 of the distribution of the on-site poten¬ 
tials. The broadened delta peak quickly comprehends the 
energy band around F = 0, reducing it to a shoulder in 
the DOS already for IF = 5. The minimum at F = — 2 is 
then filled and for IF = 10 the lowest band is already a 
shoulder only. Finally, from the three maxima only small 
shoulders remain close to the tails of the distribution. 

We then investigate the E- and IF-dependence of the 
mobility edge separating localized and delocalized states. 
The corresponding phase diagram of the CST lattice is 
depicted in Fig.|^ The critical points are mainly gathered 
from MFA and GRM. In the case of MFA we use system 
sizes from L = 20 to 60 with about 13 points for each size 
in the interval p G [pc ~ 1, Pc +1] around the critical point. 
In the case of GRM we proceed alike. We use system sizes 
from L = 3 up to 1/ = 9 and N = 4 • 10 ^. In general, 
the irrelevant scaling behavior increases with decreasing 
IF. This is more intense on the low energy side. Therefore 
critical points were calculated by GRM at higher disor¬ 
der. With ELS we only calculated the transition at fixed 
energy E = — 4 . 0 , because of the high numerical effort. 
In all cases the standard deviations of Ec and IFc are 
smaller than 0.02. As a consequence of the shape of the 
DOS, also the mobility edge is asymmetric with respect 
to the mean energy. Our results resemble the findings of 
Fazileh et al. M but improve significantly on accuracy. 
The phase diagram bears some resemblance to that of the 
fee lattice m- For low disorder strength the mobility edge 
follows the band edge. On the high energy side the mobil¬ 
ity edge is quickly drawn into the region of very high DOS 
and then it tends towards the band center with increasing 
disorder, while on the low energy side the mobility edge 
enters the band tail slowly until it approaches the local 
maximum of the DOS. Here a shoulder-like structure ap- 
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Fig. 6. Eigenfunctions of a system with 120^ UCs at fixed energy E = —4.0 for three different disorder strengths for (from left 
to right) delocalized (W = 13.5), critical (W = 14.5), and localized (W = 15.5) case. The on-site probability is visualized 
by the volume of the spheres and the indicated logarithmic color scale where A denotes the number of sites. The darkest blue 
of the color scale would correspond to a homogeneously extended state, the brightest red to a state localized on a single site. 



E 


Fig. 5. Phase diagram of the CST lattice. The area of extended 
states is enclosed by the red curve. Standard deviations for the 
data points are smaller than the red line width. The contour 
lines show the behavior of the DOS. This is divided into three 
areas (light, middle and dark grey lines; Ag = 0.004, Ag = 
0.01 and Ag = 0.1, respectively) which are separated by black 
curves at g = 0.04 and g = 0.1. Additionally, the regions of 
^ > 0.1 are shaded yellow. 


pears in the mobility edge when the local minimum of the 
DOS is traversed. The largest critical disorder strength 
Wc = 14.5 at which all states are localized is observed at 
energy E = —4.0. We choose this transition to compare 
the results of MFA, GRM and ELS with respect to the 
universal behavior by high accuracy calculations. 

For MFA we compute one state closest to energy E = 
—4.0 for each realization of disorder by diagonalization of 
the Anderson Hamiltonian 0 using a Jacobi-Davidson 
method with preconditioning |32]. Characteristic exam- 


Table 1. Number of wave functions used for each disorder 
strength depending on the number of UCs in the MFA. 


L 20-60 70, 80 90, 100 110, 120 

states 5000 3000 2000 1000 


Table 2. Number of eigenvalues used for each disorder 
strength depending on the number of UCs in the ELS. 


L 10-30 35, 40 45 50 

eigenvalues 10000 5000 4500 2500 


pies for localized, critical and extended wave functions are 
presented in Fig. For a statistical average a large num¬ 
ber of realizations for each disorder strength W is chosen 
in dependence on the system size L as given in Table 
The influence of the finite size effects can be controlled 
via the range of box sizes A used to fit the slope. Leaving 
out small box sizes reduces the irrelevant contributions to 
FSS. In general a reduction of the number of box sizes 
increases uncertainties. We use the singularity strength 
aq=o for the moment q = 0 and compare different ranges 
of box sizes between A = 0.1 and A = 0.5. The results are 
given in Table For 0.1 < A < 0.4 the irrelevant behav¬ 
ior is clearly evident and we obtain stable fits with ni > 0 
(see Fig. [^. The crossing point of the curves in Fig. 
yields the critical singularity strength aq^o{Wc) ~ 4.0. 
This value is the same as in the case of Anderson local¬ 
ization in sc lattices [lOlITT] and vibrational localization 
in fee lattices [33] . This supports the assumption that the 
multifractal spectrum at criticality depends only on the 
dimensionality and the Wigner-Dyson class [9]. 

The system-size dependent shift of the crossing point 
is nicely reproduced by the irrelevant part of the fit func¬ 
tion 0 - The histograms of the corresponding error prop- 
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Fig. 7. Behavior of aq=o for the range A G [0.1, 0.4] of box 
sizes in dependence of disorder strength W and system size. 
The lines are result of the ht with the expansion orders nr = 3, 
m = 1, TTir = 2, and mi = 0 (see second line of Table [^. 
The inset displays the scaling function and the corrected data 
according to the ht. 




Fig. 8. The histograms show the relative number of values 
(i.e. the probability distribution function) obtained for critical 
disorder strength Wc and critical relevant exponent n by the 
error propagation calculation described in Sect. [3T] applied to 
the ht result in Fig.[^ The lines show the corresponding Gaus¬ 
sian curves. 


agation in Fig. show Gaussian distributions. We there¬ 
fore assume that the ht is stable and that we can use the 
standard deviations as measure for the statistical uncer¬ 
tainties. By increasing the irrelevant expansion order n[ 
or decreasing the irrelevant behavior it is more likely to 
get deviations from Gaussian distributions, e.g. the dis¬ 
tributions are a superposition of several Gaussians. These 
results depend in particular on the initial variables of the 
minimization algorithm. For the range 0.3 < A < 0.5 the 
irrelevant behavior is signihcantly reduced. This is rein¬ 
forced by omitting all values with L < 50. These data can 
be represented by a scaling expansion without irrelevant 
exponent (ni = 0). Such regressions are very stable and do 
not depend on the initial values strongly. In general, we 
hud a higher relevant exponent u. For all results a single 
Gaussian distribution is observed. Despite the fact that 
these regressions simplify the actual behavior, they never¬ 
theless confirm the results of the fits with higher irrelevant 
influences. 

For MFA the quality Q of fit is either very close to 1 
or very close to 0 (for 0.1 < A < 0.5) indicating a problem 
in the error estimation of aq. The values of the errors can 


be overestimated in case of Q = 1 or underestimated for 

Q = 0. 

As a consistency check, we have also determined the 
probability density autocorrelation function 


C{R) = 





(29) 


with the on-site probability at the real space position 
Vi. In particular for W = 14.5 the correlation function 
shows the expected power-law behavior C{R) oc R~^ in 
dependence on the Euclidean distance R. The exponent 
c = 1.51 is close to the correlation dimension Dq =2 = 
Tq =2 = 1.39(2) obtained from Eq. Q. As in Ref. 34 the 
deviation of 8% can be attributed to unite size effects and 
the small number of samples. 

Eor the calculation with GRM we build a long bar 
with N = A • 10^ UCs. The cross section L x L is always 
square shaped with system size L ranging from 1 to 12 
UCs. We use 50 realizations of disorder for each size and 
the same disorder strengths as in the case of the MEA. 
The complex energy part r] is set to 10“^. The results 
for several expansion orders are listed in Table The 
irrelevant exponent ^ ~ 8 is higher than the value y ^ 2 
from MEA. 

The eigenenergies required for ELS are computed by 
diagonalization of the Anderson Hamiltonian Q analo¬ 
gously to determination of the eigenstates for MEA. We 
calculate the 101 eigenenergies closest to a given energy E 
for each realization of disorder. This method allows us to 
reach relatively large system sizes, it is, however, accom¬ 
panied by a high numerical cost. With increasing system 
size we decrease the number of realizations (see Table |^. 
The fit results are given in Table 

The results of the error propagation for GRM and ELS 
are very similar to those of the MEA (Eig. 1^. The results 
of all three methods are composed in Eig.^ The results 
and accuracy of MEA and ELS are very similar. GRM 
appears to be significantly more precise but agrees with 
the MEA and ELS results only in the value for the criti¬ 
cal exponent u. Eor the critical disorder strength Wc the 
distributions do not overlap indicating an extra source of 
error not captured by our statistical analysis. Therefore, 
it appears not appropriate to use the error of the mean. 
Instead, we calculate the mean and variance of the sum of 
the three distributions. In this way the variance describes 
more than the pure random part. We obtain for the tran¬ 
sition at energy Ec = —4.0 the overall critical disorder 
strength Wc = 14.474(8) and the universal critical expo¬ 
nent u = 1.565(11). While the value of Wc depends on the 
specific model the value of z/ is a universal quantity and 
agrees well with recent publications [iniin]. 

One possible source for the discrepancy of the Wc val¬ 
ues in Eig. could be that the investigated overall system 
sizes are qualitatively different: Eor MEA and ELS we con¬ 
sider large rhombohedrons (equivalent to large cubes) of 
the same extension in the three dimensions with periodic 
boundary conditions in all three directions. With GRM 
we treat quasi-one-dimensional bars of very large length 
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Fig. 9. Probability distributions of the critical parameters for 
the transition at energy E = —4.0 are described by Gaussian 
distributions. The Gaussians of the MFA are shown in Fig. 

For the data of GRM and ELS see the first line of Tableland 
the seventh line of Table , respectively. 


(which explains the high accuracy reflected in the sharp 
distribution in Fig. but relatively small cross-section, 

i.e. small extensions in the other two dimensions. That 
could lead to enhanced interference effects due to periodic 
boundary conditions in these directions and thus to the 
somewhat stronger localization so that Wc is reduced. 


5 Conclusion 

We compared the scaling behavior at the LD transition 
obtained by different methods in order to analyze the con¬ 
sistency of the obtained critical parameters and the accor¬ 
dance with universality. A sophisticated error propagation 
analysis offers the basis to estimate the accuracy of our re¬ 
sults. Generally, the differences between the methods are 
convincingly small. The averaged regular scaling exponent 
is 1.565(11) and emphasizes the universal character. This 
value is comparable with other results for the orthogonal 
universality class unHnna. However, the obtained crit¬ 
ical disorder strengths do not agree completely, but the 
deviations are smaller than the disorder-sampling inter¬ 
val. The mean value is Wc = 14.474(8). The discrepancy 
follows from unknown systematic deviations. The error es¬ 
timates describe the random deviation only. Nevertheless, 
a thorough random error analysis is important to find an 
appropriate regression of the data. 

For the phase diagram of the CST lattice we obtain 
a detailed shape of the mobility edge which is consistent 
with the results of Fazileh et al. [14] . The shape of the mo¬ 
bility edge bears some resemblance to the fee lattice m- 
We obtain a small shoulder on the lower energy side. The 
critical point at which all extended states disappear with 
increasing disorder is closer to the band center than in 
the fee lattice. The sc, the bcc and the fee lattice have 
the corresponding critical disorder strength = 16.53, 
_ 20.81, and = 26.73, respectively [TOllTS] . 

As the LD transition is (at least in this energy range) 
an interference phenomenon, it is influenced by the num¬ 
ber of nearest neighbor hopping terms N. For comparison, 
we use the relative critical disorder Wc = Wc/N. Apply¬ 
ing this to the sc, the bcc, and the fee lattice we obtain 
= 2.76, = 2.60, and = 2.23, respectively. 


This suggests that the relative value decreases with in¬ 
creasing number of nearest neighbors. With = 2.41 

the CST lattice takes a value in this range close to the 
value of the fee lattice although the number of nearest 
neighbors is the same as for the sc lattice. But as an inter¬ 
ference phenomenon, the LD transition is also influenced 
by the return probability which in turn is influenced by 
the shortest length of closed loops of nearest neighbor con¬ 
nections, which is 3 in the CST lattice and thus the same 
as in the fee lattice. 

Apart from the results, the methods have shown ad¬ 
vantages and disadvantages. In particular with GRM it is 
difficult to compute weakly disordered lattices because of 
the increasing irrelevant behavior that could be compen¬ 
sated only by higher numerical effort. On the other hand 
the measure A of GRM can be calculated directly and 
therefore the calculation is very stable. In contrast, the 
other methods which are based on eigenvalues or eigen¬ 
functions of large sparse matrices, are indirect and not as 
stable as GRM. The MFA can be used in a broad disor¬ 
der range. Considering error propagation, the analysis of 
aq by linear regression is problematic. The errors tend to 
be overestimated. Values of different box sizes may not be 
statistically independent. The ELS is simple to implement, 
but requires a large number of calculated eigenvalues or, 
to be more precise, lattice realizations. Consequently the 
system size is limited and thus restricts the accuracy of z/ 
particularly. 
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Table 3. The estimates of the critical parameters Wc, i^, and y with their standard deviations obtained from MFA for several 
expansion orders. The indices 1 and u denote the lower and upper bound of system size, respectively. Q is the quality of fit. Ny> 
and Ny are the number of points and degrees of freedom, respectively. For all transitions <^^=0 was analyzed in the disorder 
range W G [13, 5,15.5]. We set mi = 0 always. 



V 

y 

Al 

Au 

u 

Lu 

TT/j- 

Ui 

rrir 

iVp 

Ny 


Q 

14.4775(93) 

1.5694(171) 

2.1560(1466) 

0.1 

0.3 

20 

120 

3 

1 

2 

286 

274 

34.01 

1.0000 

14.4775(51) 

1.56945o1) 

2.1598(1234) 

0.1 

0.4 

20 

120 

3 

1 

2 

286 

274 

139.01 

1.0000 

14.4757(55) 

1.5653(148) 

2 . 14325252 ) 

0.1 

0.4 

20 

120 

4 

1 

3 

286 

271 

137.24 

1.0000 

14.4609(34) 

1.5889(79) 

1.8628(1000) 

0.1 

0.5 

20 

120 

3 

1 

1 

286 

275 

487.89 

0.0000 

14.4747(35) 

1.5654(135) 

2. 11725091 ) 

0.1 

0.5 

20 

120 

4 

1 

4 

286 

270 

416.87 

0.0000 

14.4705(40) 

1.5892(127) 


0.3 

0.5 

50 

120 

3 

0 

1 

208 

202 

99.41 

1.0000 

14.4825(46) 

1.5807527) 


0.3 

0.5 

50 

120 

3 

0 

2 

208 

201 

81.59 

1.0000 

14.4796(50) 

1.5906553) 


0.3 

0.5 

50 

120 

3 

0 

3 

208 

200 

80.05 

1.0000 

14.4795(51) 

1.590o558) 


0.3 

0.5 

50 

120 

3 

0 

4 

208 

199 

80.02 

1.0000 

14.46925o) 

1.5788530) 


0.3 

0.5 

50 

120 

4 

0 

1 

208 

201 

91.13 

1.0000 

14.4802(54) 

1.5784532) 


0.3 

0.5 

50 

120 

4 

0 

2 

208 

200 

80.93 

1.0000 

14.4805(53) 

1.6015558) 


0.3 

0.5 

50 

120 

4 

0 

3 

208 

199 

79.79 

1.0000 

14.4806(56) 

I. 6 O 2959 O) 


0.3 

0.5 

50 

120 

4 

0 

4 

208 

198 

79.79 

1.0000 


Table 4. Same as Table but for the A data from GRM. 


Wc 

V 

y 

U 

Lu 

Tlr 

rii 

TTlr 

Ny 

Ny 

x"" 

Q 

14.4639(8) 

1.5656(19) 

8.28(26) 

1 

12 

4 

2 

2 

300 

281 

337.09 

0.0121 

14.464o5) 

1.563459) 

8.2854) 

1 

12 

4 

2 

3 

300 

280 

336.01 

0.0121 

14.46395) 

1.56155l) 

8 . 3355 ) 

1 

12 

4 

2 

4 

300 

279 

333.11 

0.0144 

14.464o5) 

1.5627(39) 

8.28(24) 

1 

12 

5 

2 

3 

300 

277 

333.54 

0.0112 

14.46445) 

1.567457) 

8 . 2353 ) 

2 

12 

3 

1 

2 

275 

263 

328.16 

0.0039 

14.46445) 

1.5673(26) 

8.22(23) 

2 

12 

3 

1 

3 

275 

262 

328.16 

0.0034 

14.46465) 

1.566454) 

5.2o5o) 

2 

12 

3 

2 

2 

275 

259 

312.23 

0.0131 

14.46455) 

1.5676(42) 

5.18(20) 

2 

12 

3 

2 

3 

275 

258 

312.12 

0.0118 

14.46395) 

1.565759) 

8 . 2755 ) 

2 

12 

4 

1 

2 

275 

261 

320.63 

0.0069 

14.464o5) 

1.56325o) 

8 . 2955 ) 

2 

12 

4 

1 

3 

275 

260 

319.55 

0.0069 

14.464o5) 

1.5646(26) 

5.33(28) 

2 

12 

4 

2 

2 

275 

256 

307.27 

0.0154 

14.464o5o) 

1.563358) 

5.3657) 

2 

12 

4 

2 

3 

275 

255 

307.00 

0.0142 

14.46645o) 

1.5661(22) 

8.71(31) 

2 

10 

3 

1 

2 

225 

213 

245.95 

0.0603 

14.46655i) 

1.565854) 

8.7l5l) 

2 

10 

3 

1 

3 

225 

212 

245.94 

0.0549 

14.466653) 

1.566855) 

5 . 4956 ) 

2 

10 

3 

2 

3 

225 

208 

236.52 

0.0851 

14.46585i) 

1.5643(22) 

8.7455 

2 

10 

4 

1 

2 

225 

211 

238.18 

0.0964 

14.46595i) 

1.560657) 

8.7652) 

2 

10 

4 

1 

3 

225 

210 

236.86 

0.0984 

14.465752) 

1.5626(34) 

5.68(58) 

2 

10 

4 

2 

2 

225 

206 

229.53 

0.1250 

14.462i5i) 

1 . 573455 ) 


5 

12 

3 

0 

3 

200 

192 

232.48 

0.0244 

14.46185i) 

1.5689(51) 


5 

12 

4 

0 

3 

200 

191 

229.36 

0.0302 

14.46265o) 

1.568o56) 


5 

12 

3 

0 

2 

200 

193 

234.62 

0.0219 

14.46195i) 

1.566556) 


5 

12 

4 

0 

2 

200 

192 

229.68 

0.0326 


Table 5. Same as Tablebut for the 7 data from ELS. There are no irrelevant contributions taken into account, i.e. ni = 0. 


Wc 

V 

Li 

Lu 

Ur 

TTlr 

Ny 

Ny 


Q 

14.4754(31) 

1.5453(116) 

10 

50 

3 

1 

234 

228 

232.33 

0.4081 

14.4755(34) 

1.54565i8) 

10 

50 

3 

2 

234 

227 

232.30 

0.3904 

14.4761(35) 

1.5697572) 

10 

50 

3 

3 

234 

226 

227.37 

0.4620 

14.4771(36) 

1.5674(167) 

10 

50 

3 

4 

234 

225 

225.49 

0.4783 

14.4755(32) 

1.54385i6) 

10 

50 

4 

1 

234 

227 

230.54 

0.4221 

14.4825(41) 

1.5490(120) 

10 

50 

4 

2 

234 

226 

223.22 

0.5399 

14.4812(41) 

1.5614565) 

10 

50 

4 

3 

234 

225 

221.90 

0.5459 

14.4814(41) 

1.5613(166) 

10 

50 

4 

4 

234 

224 

221.83 

0.5284 
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